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Abstract Collisions of Au+Au have been studied at beam kinetic energies of 6, 
8, and 10.8 GeV/nucleon at the AGS facility at Brookhaven National 
Laboratory. Particles emitted from the collisions were momentum ana- 
lyzed and identified in a magnetic spectrometer system. Measurements 
were made at spectrometer angles in the range 14° - 59° . mt-spectra of 
protons from central collisions were analyzed to derive integrated rapid- 
ity distributions and inverse slope as a function of rapidity. The results 
are compared with a thermal model and it is concluded that there is 
either substantial transparency or longitudinal expansion at all three 
beam energies. 

Keywords: AGS, Protons, Rapidity distributions. Thermal models 

1. INTRODUCTION 

One of the expected ways of achieving the Quark Gluon Plasma phase 
of hadronic matter is to increase the matter density to 5-10 times that 
of normal nuclear matter. It has been thought that these high densities 
may be achieved in central collisions of heavy nuclei at AGS energies, 
provided that the relative motion is stopped. The question of the degree 
of stopping in head-on collisions is thus of central importance at these 
energies. Direct information pertaining to this issue may be obtained 
by studying the transverse mass spectra of protons, the majority of 
which are primordial, over a wide rapidity range to assess whether the 
observed rapidity distribution is consistent with the initial momentum of 
the projectile being converted to isotropic emission from a source at rest 
in the center-of-mass system. We find that complete stopping, defined 
in this way, is not achieved for the central Au-Au collisions at any of the 
energies studied in the E917 experiment. 

2. EXPERIMENTAL ARRANGEMENT 

The experimental arrangement is illustrated in Fig. 1.1. Beams of 
^^^Au with momenta of 6.84, 8.86 and 11.69 GeV/c per nucleon corre- 
sponding to kinetic energies of 6.0, 8.0 and 10.8 GeV per nucleon were 
obtained from the AGS at Brookhaven National Laboratory and focused 
onto a Au-target of 1 mm thickness, which corresponds to ^ 3% inter- 
action probability in the target. The trajectory of each beam particle 
was determined by a beam vertex detector consisting of four planes of 
scintillating fibers read out by position sensitive photo-multiplier tubes 
arranged in orthogonal pairs and located at 5.84 m and 1.72 m upstream 
from the target position. Further beam characterization was performed 
using beam-time-zero and halo counters also placed upstream from the 
target. Triggers for beam interactions with the target were obtained by 
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requiring that a signal of less than 75% of that expected for the full 
energy loss of a Au-beam nucleus was registered in a circular "bulls-eye" 
Cerenkov detector placed 11 m downstream from the target. 

The centrality of each beam-target interaction was derived from the 
multiplicity of particles (mostly pions) registered in a multiplicity de- 
tector array subtending a solid angle of about 6.85 sr around the target 
and/or by the total energy of the projectile remnant measured in the 
zero degree calorimeter. 

In order to determine the reaction plane orientation in peripheral 
collisions a hodoscope consisting of two orthogonal planes of 1 cm wide 
plastic scintillator slats was placed in front of the zero degree calorimeter. 
The azimuthal angle of the reaction plane determined from the average 
position of the charged projectile remnants in the hodoscope, relative 
to the beam axis will be used to study collective flow characteristics of 
peripheral collisions and the reaction plane dependence of the apparent 
source size obtained in a Hanbury-Brown Twiss analysis of pion pairs. 
Such an analysis is the subject of B. Holzman's talk at this workshop 

[!]• 

Particle spectra were obtained by momentum analysis in a movable 
magnetic spectrometer consisting of a 0.4-Tesla magnet (Henry Higgins) 
and a number of multi-wire ionization chambers used to determine the 
straight line trajectories of particles entering and exiting the magnetic 
field. In addition, a plastic scintillator wall located behind the spec- 
trometer provided particle identification by time-of-flight measurements 
relative to the Cerenkov start detector located in front of the target. 

This arrangement is capable of identifying charged pions, kaons, pro- 
tons, anti-protons and heavier nuclei. The separation of pions and kaons 




Figure 1.1 Experimental arrangement. 
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is effective up to an energy of ~1.75 GeV. Protons are separated from 
pions up to ^3.4 GeV, although there is a neghgible kaon contamination 
above -2.9 GeV. 

The acceptance of the spectrometer is also sufficient to detect corre- 
lated decay products of $ — > K^K~ ^ A pvr~, and A pir^ . The 
analysis of $-mesons and A is the subject of W.-C. Chang's talk at this 
workshop [2]. 

The present talk will concentrate on the proton spectra obtained at 
the three beam kinetic energies of 6, 8, and 10.8 GeV/nucleon, and an 
analysis to the resulting rapidity distributions and fitted inverse slopes 
obtained from the measured m^-spectra. 
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Figure 1.2 Proton mt-spectra for central (<5%) Au-Au collisions at a beam kinetic 
energy of 10.8 GeV/nucleon for different rapidity bins (laboratory system). Curves 
are Boltzmann distributions fitted to the data. Spectra for adjacent rapidity bins are 
offset by factors of ten to avoid overlapping. 
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3. RESULTS 

In Fig. 1.2, spectra of the invariant probabihty for proton emission 
per trigger event are plotted for the 10.8 GeV/nucleon beam energy as a 
function of the transverse mass m^— mg for the 5% most central collisions 
as determined from the energy deposition in the zero-degree calorimeter. 
The spectra are shown for different rapidity bins as indicated. Only 
statistical error bars are shown. Note that adjacent spectra are offset 
by factors of ten to avoid overlapping. The range in mt reflects the 
acceptance of the spectrometer in the different rapidity bins ranging from 
backwards to near mid-rapidity at ypp = 1.613. The curves represent the 
best fits to the spectra using a Boltzmann distribution, i.e. 

^ Cmtexp(-mt/r), (1.1) 



27rmt dydmt 

where C is a normalization constant and T is the inverse slope of the 
spectrum, both of which are determined from the fit to the data. We 
observe that the experimental m^-spectra are in excellent agreement with 
this shape although they could also be described almost equally well by 
a pure exponential function. 

From these fits we derive the total probability for proton emission per 
unit of rapidity, dN/dy^ which is plotted as a function of rapidity in the 
center-of-mass frame, y — ycm^ in the left panels of Fig. 1.3. The derived 
inverse slopes T are shown in the right hand panels. Data are shown for 
5% central collisions at all three beam energies, where the centrality for 
the 6 and 8 GeV/nucleon data are obtained from the the multiplicity 
array at the target position. The measured points are represented by 
solid circles, whereas reflection around mid-rapidity results in the open 
points. Error bars on the fit parameters are purely statistical and do 
not include possible systematic errors. 

We note that all three dN/ dy distributions are quite flat over the mea- 
sured rapidity range. The inverse slopes, T, show, however, a distinct 
peaking at mid-rapidity. The solid curves in Fig. 1.3 represent the ex- 
pectation for isotropic emission from a thermal source with temperature. 
To, at rest in the center-of-mass system. For such a source one expects 
that the ^/-dependence of the inverse slope is: 

T = To/coshy (1.2) 

and a dN/dy distribution of [3] 

dN ( T \ 

— oc To + 2mo — ^ + 2 — ^ exp(-mo coshy/To), (1.3) 
ay \ cosh y cosh"^ y I 
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Figure 1.3 Proton rapidity distributions in the center-of-mass system (left panels) 
and inverse slopes (right panels) are compared to a simple thermal source prediction 
for the three beam energies. The arrows indicate target and beam rapidities. 



where mo is the proton rest mass. In the right hand panels of Fig. 1.3 
we compare the rapidity dependence of the inverse slope, T, with those 
predicted by this model. The source temperature, Tq, was adjusted to 
account for the observed inverse slope at mid-rapidity. We note that this 
naive model gives a rather good representation of the observed inverse 
slopes. 

On the other hand, a comparison of the predicted distribution in 
rapidity dN/dy with the measurements (left hand panels in Fig. 1.3) 
reveals a discrepancy which clearly demonstrates that the observed pro- 
ton spectra are inconsistent with isotropic emission from a single source 
at rest in the center-of-mass system. Rather, we note that the rapidity 
distributions for protons are flat over a wide range of rapidities indi- 
cating a significant degree of either incomplete stopping or longitudinal 
expansion at all three beam energies. 
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Of course, the naive model shown here also disregards the possible 
effects of radial expansion in the fireball. It has been shown[3], however, 
that, within a wide range of parameters, there is a strong anti-correlation 
between the radial expansion velocity, 'Uq, and the source temperature, 
Tsource^ such that it is impossible to disentangle their relative values from 
fits to spectra of a single particle species e.g. protons. In the present 
analysis we have therefore chosen to use only a single parameter, namely 
the apparent source temperature, Tq, keeping in mind that its value does 
not necessarily represent the true temperature of the source formed in a 
central Au-Au collision. 
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Figure 1.1^. Proton rapidity distributions (top) and inverse slopes (bottom) for Au-Au 
collisions at 10.8 GeV/nucleon beam kinetic energy compared with previously pub- 
lished results from the E866 and E877 experiments at the AGS. The arrows indicate 
target and beam rapidities. 

In Fig. 1.4 we compare the dN/dy (top panel) and inverse slopes T 
(bottom panel) from our experiment to results from the E866 [4] and 
E877[5] experiments at the AGS. At mid-rapidity there is good agree- 
ment between the present data for the dN/dy distribution and that from 
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the E866 experiment, but a relatively small discrepancy between the two 
data sets is apparent at less central rapidities. The source of this dis- 
crepancy is presently being investigated. The E877 data were measured 
in the extreme forward rapidity region but overlap with the present (re- 
flected) data in the rapidity region y — ycm = 0.8 - 1.1. Here, there 
appear to be substantial discrepancies between the data sets. The gen- 
eral trend of the three data sets is, however, clear. There is a wide range 
-0.7< y — ycm < 0.7 around mid-rapidity, where the dn/dy distribution 
is essentially flat, followed by a monotonic decrease on either side. We 
note that this observed range of essentially constant dN/dy is inconsis- 
tent with thermal emission from a stopped source (solid curve), as well 
as early predictions of the Relativistic Quantum Molecular Dynamics 
model[6] (RQMD vl.08, [5]) (solid histogram), both of which are too 
sharply peaked at mid-rapidity. 

The slopes, however, appear to be in reasonable agreement with the 
1 / cosh y dependence expected from the thermal model (solid curve) al- 
though this may be fortuitous since the rapidity distribution clearly show 
that there is either a significant amount of transparency or longitudinal 
expansion at these energies which violates the assumption of a thermal 
source at rest in the center-of-mass system. 

4. CONCLUSION 

An analysis of m^-spectra for protons emitted in central Au-Au colli- 
sions at beam kinetic energies of 6, 8 and 10.8 GeV/nucleon in terms of 
Boltzmann distributions has been carried out, and the resulting dN/dy- 
distributions and inverse slopes derived. They are compared to a simple 
thermal model assuming isotropic emission from a source at rest in the 
center-of-mass system corresponding to complete stopping of the col- 
liding Au nuclei in central collisions. We find that the rapidity distri- 
butions dN/dy are substantially wider and essentially constant around 
mid-rapidity although the inverse slopes exhibit the expected 1 / cosh y 
dependence on rapidity. We interpret this as a manifestation of incom- 
plete stopping or longitudinal expansion of the entrance channel mo- 
menta at all three beam energies. 
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